NONINTERSECTING BROWNIAN MOTIONS ON THE HALF-LINE AND 
DISCRETE GAUSSIAN ORTHOGONAL POLYNOMIALS 
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Abstract. We study the distribution of the msiximal height of the outermost path in the model of 
TV nonintersecting Brownian motions on the half-hne as A'^ ^ oo, showing that it converges in the 
proper scahng to the Tracy- Widom distribution for the largest eigenvalue of the Gaussian orthogonal 
ensemble. This is as expected from the viewpoint that the maximal height of the outermost path 
converges to the maximum of the Airyj process minus a parabola. Our proof is based on Riemann- 
Hilbert analysis of a system of discrete orthogonal polynomials with a Gaussian weight in the double 
scaling limit as this system approaches saturation. We consequently compute the asymptotics of 
the free energy and the reproducing kernel of the corresponding discrete orthogonal polynomial 
ensemble in the critical scaling in which the density of particles approaches saturation. Both of 
these results can be viewed as dual to the case in which the mean density of eigenvalues in a random 
matrix model is vanishing at one point. 



1. Introduction 

1.1. Nonintersecting Brownian motions on the half-line. Consider a model of N noninter- 
secting Brownian motions {bj{t)}jLi which remain non-negative for < t < 1 and whose initial 
and terminal points are at zero. That is, 

6i(0) = 6i(l) = 62(0) = 62(1) = ••• = &^(0) = &w(l) = , 

<bi{t) <b2{t) < ■■■ <bN{t) for 0<t<l. ^ ' 

There are two standard ways to enforce the condition that the Brownian motions remain non- 
negative: an absorbing wall and a reflecting wall at zero. The transition probability for a single 
Brownian motion with an absorbing wall at zero to pass from y to x over the time interval t is 
given by 



Pabs{t,x\y) 



1 



exp 



exp 



(y + xY 



(1.2) 



2t J ^ V 2t 

and the transition probability for a single Brownian motion with a reflecting wall at zero to pass 
from y to X over the time interval t is given by 



Pref{t,x\y) 



1 



exp 



2t 



+ exp 



2t 



(1.3) 



For the positions at time t £ (0, 1) of the nonintersecting Brownian motions with an absorbing 
wall at zero, we will use the notation 



< 6' 



(^^)(t)<6f^)(t)< 



(1.4) 



'1 ^ "2 

The superscript BE stands for Brownian excursion, which is the name for a Brownian motion with 
an absorbing wall which is conditioned to return to its starting point. For the positions at time 
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t G (0, 1) of the N nonintersecting Brownian motions with a reflecting wall at zero, we use the 
notation 



0<b[''\t)<bi''\t)<---<b^^\t), 



(1.5) 



where the superscript R stands for reflecting. The ensembles of nonintersecting Brownian motions 
can be derived from the transition probabihties (1.2) and (1.3) using the Karlin- McGregor formula 
[30]. Even though it seems Hke a degenerate condition to force all of the Brownian motions to 
begin and end at zero, it is possible to define these models of nonintersecting Brownian motions 
with such a condition by starting with a model for which the starting and ending points are all 
distinct and positive, and taking a limit as they go to zero. Let us now give a precise definition of 
the two models of nonintersecting Brownian motions in terms of their transition probabilities. See 
[34] for a derivation of these transition probabilities in the absorbing case. The reflecting case is 
similar. Introduce the notations 



Xfc — {xk,l,Xk,2, ■ ■ ■ , Xk,N), — (a^fc,l) a^fc,25 • • • > ^k,N)^ 

and let 

^(x)= n (xk-xj), 

l<j<k<N 

be the Vandermonde determinant. Define the probability density functions 

(t(l-t))-^'-^/'2^/2 



(1.6) 
(1.7) 



N 



{BE) 

Po 



(«,xi) 



^^/^nf=o'(2j + l)! 



(A(x?))^ n4. h-p 




(t(l-t))-"^^"/V/^ , / 1 f , ^ 



(1.8) 



^^/^n;=o(2i)! 



and 



p(^^)(tl,Xi;t2,X2) = 



l-t2 

X det 



N 



A(x|) 
A(x?) 



exp 



Xl,j 



i-t2 i-h 



Pabsit2 - tl,X2,j\xi k) 



;j(^)(il,Xi;t2,X2) 



l-tA^^-^/^A(xi) ^ 



N 

j,k=l 

N 



2 ^ V 1 - t2 1 - 



X det 



Prefih - ti,X2,j\xi^k) 



N 

j,k=l 



(1.9) 

Let us adopt the convention that bj{t) with no superscript refers to either the model with the 
reflecting or absorbing wall at zero. For any time t G (0, 1), the ordered particles (6i(t), . . . , b]sf{t)) 
must lie in the region in 

W;v = {{bi, . . . , 6„) : < 6i < 62 < • • • < b^v}- (I-IO) 
The meaning of the probability density functions above is the following. For some sequence of times 

<ti < ■■■ <tK <l, (1.11) 
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and some sequence of regions cWn, k = 1, . . . ,K we have 
{hitk),b2{tk),...,bNitk)) G Afc, k = l,...,K 



/ ••• / Po(ii,xi) JJp(tfc,Xfc;ifc+i,Xfc+i) (ixi---dxx, 
^Ai JAk Vfe=l / 



(1.12) 



where Xjt is the vector of integration variables corresponding to the region A^, po = P(^^^ (resp. 

Pq^^), and p = p^^^^ (resp. p^^^) for the model of nonintersecting Brownian motions with an 
absorbing (resp. reflecting) wall at zero. 

From (1.8) it is immediate that for fixed t € (0,1), the particles b- (t) are distributed as 
the eigenvalues of a random matrix from the Laguerre unitary ensemble (see e.g., [24]). As such, 
the largest particle at each fixed time, in the proper scaling limit, is distributed according to the 
distribution which describes the largest eigenvalue in the Gaussian unitary ensemble (GUE) of 
random matrices as the size of the matrix tends to infinity. In this paper we study the distribution 
of the random variable 

max 6jv(i), (1-13) 

0<t<l 

the maximal height of the outermost path for both the absorbing and reflecting case. To obtain a 
formula for the distribution of this random variable, one can apply the Karlin-McGregor formula 
in the afiine Weyl alcove of height M 

= {(6i, . . . , 6„) : < 6i < 62 < • • • < biv < M}. 



(1.14) 



In the absorbing wall case, the formula obtained is 



max b^^^\t) < M 



\0<t<l 



2-N/2^2N'^+N/2 



N 



E (A(x^))Mn-. 



-vr 



exp 



TV 



(1.15) 



where 

X = {X1,X2, . . .,Xn)- 

Notice that this is the Hankel determinant 



2-N/2^2N'^+N/2 



det 



xe 



-7r2/2M2, 



N 



j,k=l 



(1.16) 



(1.17) 



This formula first appeared in the paper [39] , in which the authors use a path integral technique to 
derive it. A derivation using the Karlin-McGregor formula appeared soon after in [34]. See also [22] 
in which an equivalent formula is derived from lattice paths. In the case that there is a reflecting 
wall at zero, the formula is 



2-Ar/2^2Ar2-3Ar/2 
MJV(2^-i)7V!n£o^(2A;)! 



(2 ^ 



Xj 



xe{Z-i/2}' 



(1.18) 

To our knowledge, the formula (1.18) has not appeared in the literature before, but it can be derived 
in a manner similar to those used in [39] and [34] to derive (1.15). 
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In the case of an absorbing wall, this model was first introduced in the papers [3f] and [26], and 
is often called the model of nonintersecting Brownian excursions, or noncolliding Bessel bridges. It 
is also sometimes referred to as "watermelons with a wall," although this sometimes refers to the 
discrete time and space (simple random walk) version as well [26] . See [26] and [32] for a derivation 
of this model as a scaling limit of an ensemble of simple random walks conditioned not to intersect 
and to stay positive. 

Our analysis of (1.15) and (1.18) is based on analysis of the discrete Gaussian orthogonal poly- 
nomials {Pk°'\x)}'^^Q and their normalizing constants {/i^"'*}^q defined via the orthogonality con- 
dition 

^ Pt\x)PtHx)w{x) = hl^^Smn, w{x) = exp i , Pt\x) = x^ + ■ ■ ■ . (1.19) 

A routine calculation (see [13]) shows that (1.15) and (1.18) can be written as 



{ max b^^^\t) < m) = ^ \\ , n 4ti 

/ , , \ r,-N/2^2N'2-3N/2 ^"1 , , , 

P( max 6if)(t) < M ) = \ \ ^ / TT 4f \ 



(1.20) 



In a recent paper of Forrester, Majumdar, and Schehr [25], an analogy beween nonintersecting 
Brownian excursions and Yang-Mills theory on the sphere is made, and the authors use some non- 
rigorous methods from gauge theory ([12], [20]) to deduce that the maximal height of the outermost 
path in this ensemble is, in the proper scaling limit, distributed as the largest eigenvalue in the 
Gaussian orthogonal ensemble (GOE) of random matrices. A main result of this paper is a rigorous 
verification of this fact. In order to state this theorem, let us review the Tracy- Widom distributions 
which describe the location of the largest eigenvalue in the classical random matrix ensembles. 
These distributions may be described in terms of the Hastings-McLeod solution to the Painleve II 
equation. The homogeneous Painleve II equation is the second order nonlinear ODE 

q"{s) = sq{s) + 2qisf. (1.21) 

The Hastings-McLeod solution to this equation [27] is characterized by its behavior at positive 
infinity. In particular, it is the solution satisfying 

q{s)=Ai{s){l + o{l)) = ^-^-^{l + o{l)j as s ^ +oo, (1.22) 

where Ai is the Airy function. The distribution functions Ti and T2 are defined as 

Ti{x) = F{x)E{x) , T2{x) = F{xf , (1.23) 

where 

E{x) = exp q{s)d^ , F{x) = exp R{s)d^ , R{x) = ^ q{sfds, (1.24) 

and q{s) is the Hastings-McLeod solution to the Painleve II equation. The function F\ describes 
the distribution of the largest (or smallest) rescaled eigenvalue in GOE, while J-2 describes the 
distribution of the largest (or smallest) rescaled eigenvalue in GUE (see [41], [42], [43]). We now 
state our main theorem. 
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Theorem 1.1. (Distribution of the maximal height of the outermost particle) Consider either of 
the models of nonintersecting Brownian motions described in (1-8) and (1-9). The maximal height 
of the outermost particle has the limiting distribution 



lim P 



211/6^1/6 [ ^ _ ,J2n\ ^ 



= J'i{k), (1.25) 



where F\, defined in (1.23) and (1.24), is the limiting distribution function for the location of the 
largest eigenvalue in the Gaussian orthogonal ensemble of random matrices, and bN{t) is either 
6Sf^)(t) orb^^\t). 

This theorem is widely expected from the point of view that the distribution of the uppermost 
curve in the model of N nonintersecting Brownian bridges should converge (after rescaling and 
recentering) to the Airy2 process, which was first introduced in [36]. In the case of an absorbing 
wall, the framework to prove this convergence at the level of finite dimensional distributions was 
given by Tracy and Widom in [44], although they stopped just short of stating it as a theorem 
(their main interest in that paper was the asymptotics of the bottom curve). It is known that 
the maximum of the Airy2 process over a continuum of times is given by the Tracy- Widom GOE 
distribution. Such a result was first proved by Johansson [29] by first proving a functional limit 
theorem for the convergence of the polynuclcar growth (PNG) model to the Airy2 process and 
using connections between PNG and the longest increasing subsequence of a random permutation 
found by Baik and Rains [4]. A more direct proof was recently given by Corwin, Quastel and 
Remenik [11]. See also [35]. Thus Theorem 1.1 could be proved by establishing the functional 
convergence of the top curve to the Airy2 process. In fact, for the absorbing boundary case, given 
the finite dimensional convergence implied by [44] the functional convergence follows from a result 
of Corwin and Hammond [10], who showed that finite dimensional convergence of a line ensemble 
implies functional convergence given some fairly mild local condition. In this sense, at least in the 
absorbing case, Theorem 1.1 is not new, but here we give an alternative direct proof by analyzing 
the formula (1-15) asymptotically. Moreover, the analysis is uniform for both the absorbing and 
reflecting boundaries. The rigorous result for the reflecting boundary case does seem to be new. Our 
proof is based on the asymptotic evaluation of the formulas (1.20) by Riemann-Hilbert methods. 

Let us note here that in the paper [25] the authors give expressions similar to (1.15) for the 
normalized reunion probabilities for nonintersecting Brownian motions with periodic boundary 
conditions and with reflecting boundary conditions, which are shown to correspond to partition 
functions of 2-d Yang- Mills theory on the sphere with the gauge groups U(A'') and S0(2A'^), respec- 
tively. These expressions do not have a probabilistic interpretation, but can also be expressed in 
terms of discrete Gaussian orthogonal polynomials and their asymptotic evaluation is a straight- 
forward application of Theorems 1.1 and 1.3. 

1.2. Discrete Gaussian orthogonal polynomials. For asymptotic analysis, it is convenient to 
use a rescaling of the polynomials (1.19). Consider the infinite regular lattice of mesh 1/n, 

1 r 

2' 2 

and the polynomials orthogonal with respect to a discrete Gaussian weight on this lattice. More 
specifically, consider the system of monic polynomials {P^j {x)}'^q and the normalizing constants 
{/i|j"-}^o satisfying the orthogonality condition 



(1.26) 



T 
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a>l 



Figure 1. The equilibrium measure density for discrete Gaussian orthogonal poly- 
nomials in the subcritical case a < 1, the supercritical case a > 1, and the critical 
case a = 1. 

As usual, Pn^'ki^) is a polynomial of degree k, and a > is a positive parameter. As the mesh of the 
lattice goes to zero these polynomials converge to the (monic and rescaled) Hermite polynomials. 
The polynomials P^^{x) and the normalizing constants h^j^ depend on the parameter a. To 
highlight that dependence, let us write 

hit - h^:lia). (1.28) 

The relation to the polynomials (1-19) is 

Pi")(nx) = n'pi^iix) , = n"^^'h^:l (a), M = ^. (1.29) 

The main distinguishing feature between the asymptotic analysis of discrete orthogonal poly- 
nomials and that of continuous ones is the phenomenon of saturation. If we denote by /i„ the 
normalized counting measure on the zero set of the polynomials Pn"n{x), it is known that, as 
n — )• oo, /in converges to a probability measure with finite support and piecewise smooth density, 
known as the equilibrium measure. It is a general fact that for any system of polynomials orthog- 
onal with respect to a measure which lies on a discrete subset of M, call it D, all zeroes of the 
polynomials are real and there can be no more than one zero between two consecutive nodes of D. 
This leads to an upper constraint on the distribution of zeroes. For polynomials orthogonal with 
respect to a continuous weight, there is no such constraint. 

If the discrete orthogonal polynomials are such that the equilibrium measure does not approach 
this upper constraint, then their asymptotic properties match those of a corresponding continuous 
system. If the upper constraint is active, then they do not, see [3], [7]. In the case of the discrete 
Gaussian orthogonal polynomials described in (1.27), the mesh of the lattice is 1/n and thus the 
the upper constraint on the equilibrium measure is that it should have a density which is no greater 
than 1. In the case that the parameter a is greater than 1, this upper constraint is realized, meaning 
that there is an interval on which the density of the equilibrium measure is identically 1 , see Figure 
1. 
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In [3] and [7] the authors present the asymptotic properties of very general classes of discrete 
orthogonal polynomials on the real line assuming some "regularity" condition on the equilibrium 
measure. The proof of Theorem 1.1 requires that we explore the critical case in which the upper 
constraint is approached in a double scaling limit, which is not considered in [3] or [7]. For the 
continuous weight case, a similar double scaling limit was studied in the seminal paper of Baik, Deift, 
and Johannson [1] in the context of the longest increasing subsequence of a random permutation. 
In the context of random matrix theory, this type of double scaling limit appears when the limiting 
mean density of eigenvalues in a random matrix model vanishes at one point, see [5], [9]. In the 
present paper we adapt the analysis to a discrete weight where the lower constraint is replaced by 
the upper constraint. A similar double scaling limit was considered recently by Baik and Jenkins 
[2] for a system of discrete orthogonal polynomials on the circle when the upper constraint is about 
to be active. 

The orthogonal polynomials (1.28) satisfy the recurrence relation (see e.g., [40]) 

In the case that a = or a = 1/2, the lattice is symmetric about zero and A^^j. is zero. The 
polynomials P^^j of course depend on the parameter a, and if we show that dependence by writing 
Plf-{x;a), then we have the lattice rescaling relations 

(1.31) 

where 

e± = 1 ± -. (1.32) 

n 

A basic physical model described by these orthogonal polynomials is the discrete orthogonal 
polynomial ensemble with Gaussian weight, which is a discrete version of GUE. This ensemble is 
described as the probability distribution on ^.-tuples of points A — (Ai, A2j • • • , A^) £ (-^^n,a)" 

P(there are particles at each of the points Ai, . . . , A^) = (z^^'^^'^n JJ(Aj - A^)^ JJ 



k>j 3=1 



where 



(1.33) 



" -nV{\j) 



and 



Ae(L„,a)" l<j<fe<n j=l fc=0 

V{x) = ^x'. (1.35) 
As proved in Appendix C, the partition function Z^^^^^ satisfies the deformation equation 

£ log Zi^OPE) ^ ^!|! j ' ^(oO (5(cO_^ + ^H^^ + (^(a) + . (i.gg) 

The deformation (1.36) is one of the isospectral flows in a general system known as the Toda 
lattice hierarchy, see [15], [17] and references therein. For a derivation of the Toda lattice from 
the Ricmann-Hilbert problem for orthogonal polynomials, see [3], and for a broad description of 
differential equations related to orthogonal polynomials, see [21]. Let us note that there also 
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exist similar deformations with respect to other parameters in the weight which yield closed form 
expressions for the first logarithmic derivative of the partition function, see [14], [28], [33]. 
If we define the free energy as 



(1.37) 



then (1.36) reads as 

do? 



\K,n-2ia) K^l^_i{a) /i;,,^_i(a) J 



(1.38) 



In light of the deformation equations (1.36) and (1.38), let us write 

ZiDOPE) ^ ZI^OPE)^^^ ^ piDOPE) ^ Fi^OPE)^^^_ (^ gg) 

All correlation functions for this ensemble can be written in terms of a reproducing kernel which 
is defined in terms of orthogonal polynomials. Introduce the V'-functions 



(a) 
n.k 



-ii\'(,i-)/2 



and the Christoffel-Darboux kernel 



n-l 



k=0 



(1.40) 



(1.41) 



(1.42) 



Then the m-point correlation function iim „(Ai, . . . A„) is defined by the formula 

Rm,n{^l, ■ ■ ■ Am) = det (K„(Afc, A/))^^^-^. 

The proof of Theorem 1 . 1 requires an asymptotic formula for (a) in the scaling limit a — >■ 1 
as n — > cxD. In order to state that expansion, let us first fix some notations. Let the parameter s 
be defined in terms of a as 

2/3 

s = s{a; n) = — 



(37rn) ^zi - p{v)dv^ 



where 



zi = —Va-1, p{n) = ^J-^-'n 



ira 



TT^a 



One may check that as a — t- 1, 

. = 2V3„2/3 1^(1 _ + _ ^)2 + _ a)^)) 



and thus if 
then 



a = 1 - xn"^ + 0(n~^) , < 5 < e, 



lim s 
n->oo 



( (sgna;)oo if < 5 < 2/3 
2^/^x if J = 2/3 
if 2/3 < 6. 



(1.43) 

(1.44) 

(1.45) 
(1.46) 

(1.47) 
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Let q{s) be the Hastiiigs-McLeod solution to the Painleve II equation, and R{s) be as defined in 
(1.24). We then have the following proposition. 

Proposition 1.2. (Asymptotics in discrete Gaussian orthogonal polynomials) Consider the or- 
thogonal polynomials defined in (1.27) such that (1 — a)n^/^ remains bounded as n ^ oo. Let s be 
defined in terms of a as in (1.43). The normalizing constants for the orthogonal polynomials (1.27) 
satisfy, as n oo, 

9 / 1 \ " / 9^/3 9I/3 \ 



, X 1 / 92/3 9I/3 \ 



(1.48) 



where 

T„(s) = R{s) - (-l)"cos(27ra)g(.s) 



(1.49) 



Un{s) = R^is) - (-l)"cos(27ra)(g'(s) + 2g(s)i?(s)) - ^^(s) sin2(27ra). 
The recurrence coefficients A^^j^_-^ satisfy 

<L» = (^'^(») + ^ - ''c-'')) ■ (-0) 

Let us also note a formula for ratios of normalizing constants which will be useful in the proofs 
of Theorems 1.1 and 1.3. Introduce the notations 

s+ = s(a^+; n + 1) , s_ = s(a^_; n — 1) , s = s{a;n), (l-^l) 

where s{a;n) is defined in (1.43), and ^± is as defined in (1.32). Assume that a = 1 — xn"^/^ for 
X G M. A direct application of Proposition 1.2 gives 



2I/3 



(?7„(s_) + C/„(s) - 2T„(s)T„(s_)) + 0( 



-1^ 



Prom (1.45) we see that 



2^/3 , \ 

+-^(C/„+i(s+) + [/n+i(5) - 2T„+i(s)T„+i(s+)) + 0(n-^) 1 . 

(1.52) 



22/3 

s - s± = ±— + O(n-i) , (1.53) 

and thus we can rewrite (1.52) as 

h^^lia) 1 / , 2V3 



h^:l,.+i{^+a) ^ ^_ , 2^/3 



1 + ::278 (^n(s) + C^n(s) - Tr,{sf ) + 0(n-^) j , 

(1.54) 
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In fact, it easy to see, using R'{s) = —q^{s), that Un{s) — Tn{s)'^ = T^{s), and thus we have 

,o) , .„ . . (1-55) 

'>»+l,„+l(f+») 



A study of the asymptotic properties of the discrete Gaussian orthogonal polynomials in the 
critical scaling above naturally leads to asymptotic results for the discrete Gaussian orthogonal 
polynomial ensemble (1.33) in the critical scaling such that the distribution of particles is ap- 
proaching saturation. The theorems below give such results, and we note that they are nearly 
identical to the results in [5], [6], and [9], which concern a random matrix model for which the 
distribution of eigenvalues is vanishing at a single point. 

We first compare the free energy in the discrete Gaussian orthogonal polynomial ensemble with 
that of the Gaussian unitary ensemble. The free energy of GUE is defined as 



l<j<k<n \ 3=1 J 



dxi ■ ■ ■ dXn- 

(1.56) 



+ 0{n-^), (1.57) 



Theorem 1.3. (Free energy in the discrete Gaussian orthogonal polynomial ensemble) Let the 
parameter a be such that (1 — a)n^/^ remains bounded as n — )■ oo. The free energy of the discrete 
orthogonal polynomial ensemble, F^^^^\a), defined in (1.34) o-'nd (1.37), satisfies as n oo 

Fi^OPE)^^^ _ piGUE) = _ 1 log _ J_ log k(2V3^2/3(i _ 

where J^2{x) is the Tracy-Widom distribution function associated with the largest eigenvalue of 
GUE, and2<5< 7/3. 

We would now like to describe the Christoffel-Darboux kernel (1.41) close to the origin as n ^ oo. 
In order to state the theorem, we need to fix some notation. Let $^(C; s) and ^^(C; *) be defined 
via the system of differential equations 

5_ (C; s)\ / ACq AC + S + 2g2 + 2r\ (^^ (C; s) 

dC \<^'{C; s) J [-Ae -s-2q' + 2r -ACq J {^HCl s), 

d_f^\C,s)\_fq C\f<i>HC;sr 

ds \^\Q:s)) - \-C -q) \<A>\Q-.s)^ 
satisfying the properties that <I>^(C; s) and <I*^(C; s) are real for real C, and s, 

$i(-C;s) = $^(C;s), $2(-C;s) = -$2(C;s), (1.59) 
and have the real asymptotics 



(1.58) 



$nC; s) = COS + + 0(C-'), «>'(C; s) = - sin {^C" + s^j + o{Q-^) 



(1.60) 



as C — ^ ±oo. These are the so-called psi-functions associated with the Painleve II equation. We 
then define the functions 

^^{C,s) = ^\C,s) + i^\C,s), ^2{C,s) = ^\C:s)-i^\C,s), (1.61) 
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and the critical kernel 



.cTitt „^ _ -*i(^; s)^2{v] s) + $2(u; s)$i(v; s) 



Let us also give two other expressions for K'^^^ which may be useful for analysis (see e.g., [9]), 



(1.63) 



As shown in [5] and [9], this is the limiting correlation kernel of a random matrix model in the 
case that the limiting distribution of eigenvalues vanishes at a single point. We have the following 
expression for the kernel (1.41) as n — >■ oo. 

Theorem 1.4. (Limiting Christoffel-Darboux kernel near the point of saturation) Let Kn{x,y) be 
the Christoffel-Darboux kernel defined in (1-41) for the discrete orthogonal polynomial ensemble 
(1.33) with potential V{x) given by (1.35). Consider the scaling 

1 r -2/3 kn-a u run-a v 5,3 
a = l-Ln x = ~ — y= ~ ~r7s' c = 7r2 (1.64) 

where kn and rUn are integers. Then, for u ^ v, 

„2/3 

lim {-\fn+mn + l'J_K^^^^y) = K'^^\u, V' s^) , (1.65) 
n-^oo c 



where 

The diagonal terms satisfy 



= 2^/^L. (1.66) 



^2/3 



lim 1 - Kn{x, x) = K'"\u, u; s^), (1.67) 

n-^oo C \ J 

where K'^^*{u,u; s^o) is obtained from (1.62) using L'Hospital's rule or directly from (1.63). 

Theorem 1.4 may have some application to nonintersecting Brownian excursions. In a recent 
paper of Rambeau and Schehr [37], the authors derive a formula for the joint distribution of the 
maximal height of the outermost path and the time at which it occurs. Their formula can be 
written in terms of the Christoffel-Darboux kernel (1.41) [38], and thus Theorem 1.4 may be of use 
in the asymptotic analysis of this joint distribution. In fact, very recently Schehr gave a limiting 
formula for this joint distribution in the critical scaling which involves the Painlcvc II psi- function 
[38]. The argument of [38] is based on a differential Ansatz, and it would be interesting to see if 
one could give a rigorous verification of that result using Theorem 1.4. 

The rest of this paper is organized as follows. In section 2, we derive an integral formula for 
the distribution of the random variable maxo<f<i 6Ar(t) in terms of discrete Gaussian orthogonal 
polynomials and use Proposition 1.2 to evaluate it in the large limit, which proves Theorem 
1.1. In section 3, we prove Theorem 1.3 in a similar way. In section 4, we present the steepest 
descent analysis of a Riemann-Hilbert problem for the discrete Gaussian orthogonal polynomals by 
the method of Deift and Zhou [18] , and in section 5 we explicity compute the first two error terms 
of this analysis in the critical scaling limit. Finally, in section 6 we use the results of sections 4 and 
5 to prove Proposition 1.2 and Theorem 1.4. 
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2. Proof of Theorem 1.1 

2.1. Integral formula for the distribution of the maximal height of bN{t). We would 
like to study the double scaling limit of (1.20) as iV — )■ oo, and M = ■\/2N + kN~^/^ for some 

k e M.. With that in mind, we scale M as M = and will study the limit as iV ^ oo and 

a = 1 - LN-^/^ + 0{N-^/^). Then formulas (1.20) become 

n ^2fe+i' 

fc=o ^2.1) 



max b^^^\t) < M 
\o<t<i •'^ 



2iV2^Jv2+iv/2 Yl^-^\2k + 1)! 

^N'^-N/2^2N'^-ZN/2 



max 6\7^(i) < M , , , , i 

Vo<t<i ^ ^ ' ) 2^'-^iV^'-^/2n£o (2^)! fe=o 



,(1/2) 



As proved in Appendix C, the products of normalizing constants in (2.1) satisfy the deformation 
equations 



\ A:=0 / 



N-l 



log n h 



(1/2) 
2fc 



fc=0 



It follows that, if we denote 
Fif'"*(o)=log 



o-?s^r'(*)<v/f) 



Ff)(a) = log 



7r2^ 


\ "'2Ar+l 




"'2N-1 


7r2^ 


2 ^(1/2) 


4Af, 


1 Am 

"■2N-2 



(2.2) 



P ( max b^^\t) <\l — 



then we have the deformation equations 

2Ar2 + iV / ^2 ^2^(0) 



do? 



2a2 



+ 



vr 
4iV 



2Ar+l 



(0) 

2iV-l 



0<t<l 



2N^ -N 



2N 



2.3) 



2a2 



+ 



iiV 



2 \ 2 , (1/2) 



2N 



h. 



In the scaling of M described above, we can write the orthogonality condition (1.19) as 

^2 



(1/2) • 
■27V-2 

(2.4) 



^ Pf)(nx)Pir)(nx)e-^W = 4°)5fc^, V{x) 



,.5) 



This is the same orthogonality condition as (1.27) with a a^, and thus using (1.31) we can write 
the formulas (2.4) as 



d2 



» p(SE), X _ 2N'^ + N ^ _4^^^2HN+3 ^2N+l,2N+M^+) 
'J^^N W — KZ9 I- TT iV 5^ 



da2 



^(R) 
iV 



(a) 



2a2 



2N^ -N 



h 



fa) 

2Af,2Af-lV"J 



2a2 



+ ^4^2^-4iV+3. 



,(1/2) 
hN,2N 



(a) 



(1/2) 

'2Ar-l,2iV-2 



-2K-) 



1 
2N 



(2.6) 



Notice that as a — 0, M — 00. Since M is typically close to V2N, it is reasonable to assume that 
both Fj^ (a) and (a) go to zero very quickly as a ^ 0. One might easily guess the following 
lemma. 
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Lemma 2.1. 



Furthermore, as N oo, 



.,..^r(«) 

a=0 



= 0. (2.7) 

a=0 



'2 



Fr\a) 



= 0{N-^), ^Fif^(a) =0{N-^). (2i 



a=Q 

We leave the proof of this lemma the Appendix A. From (2.7) and (2.6), we easily obtain the 
following integral representations for F^^\a) and F^\a). 

Proposition 2.2. The functions F^^\a) and F^\a), defined in (2.3), have the integral repre- 
sentations 

pa ru pa pu 

Fjf'^>{a)= / / G^^^\r)drdu, Fj;^> (a) = / / G^^\r)drdu, (2.9) 
Jo Jo Jo Jo 

where 

r(BE)(^^ _ 2Ar2 + Ar 4 2^4iV+3 ^2iV+l,2iV+lK+) 

h2N,2N Ar) ^^^^^^ 

^ ' 2r2 ^" i,{i/2) f t \' 

"'2Ar-l,2Ar-2V s-j 

2.2. Evaluation of the integrals (2.9). We would like to evaluate the integrals (2.9) in the limit 
as iV ^ 00, and 

a = 1 - + 0(A^-^/2). (2.11) 

Let us write (2.9) as 

Fif'^Ha) = + + if F(fV) = 4"^^ + /i^^ + (2.12) 



where 



^) = /'""^ ' rG^^^\r)drdu, if ^) = /'""^ ^ G(^^)(r)drdn, 

Jo Jo Jl-N-^ Jo 

[•a i-u rl—N~^ pu 

/f^)=/ / G(^^)(r)drd^, /f)= / / G(«)(r)drdn, (2.13) 

Jl-N-^ Jl-N-^ Jo Jo 

pa pl—N~^ pa pu 

l[^> = / G^^\r)drdu, 4 = / / G^^\r)drdu. 

Jl-N-^ Jo Jl-N-^ Jl~N-e 

for some < e < 2/3. 

Consider first Iq. We need a large n formula for the normalizing constants h^^{r) when r < 

1 — n~^. In this case the asymptotics of h'^'^j,{r) match the asymptotics of the corresponding system 
of continuous orthogonal polynomials, the (monic and rescaled) Hermite polynomials. We have the 
following lemma, whose proof is given in Appendix B. 

Lemma 2.3. Let h^j.{a) be defined as in (1.27) and (1.28). Let the parameter a he such that 
a < 1 — for some < £ < 2/3. Then as n 00, 



H / ^-l _ 1 ( 2 \ri ( . 1 , 1 , 139 _4, 



(2.14) 
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Combining these, we find that, as A/" — >■ oo. 



"'2N,2N-1V ) ^ 



(2.15) 



(2.16) 



;,(V2) /„t ^1 7r4r2e2 V 2N 24iV2 i2Af3 

"'2N-l,2N\' 'i-) ^ 

Inserting this asymptotic formula into the integrand of l'^^^ and /q^^ gives 
Similarly, we find that 

Ip^) = 0(iV-2-^) , if^ = 0{N-'^-') . (2.17) 

(BF) ( R) 

Wc arc left to evaluate I2 and • These integrals are in the regime in which Proposition 
1.2 is valid. Let us write r = 1 — xN~'^^^, so that as r varies from 1 — A'^"'' to a, x varies from 
]\j2/3-e ^ Applying equation (1.55) with n = 2N, we obtain 



^2JV+l,2Af+l(^^+) 
"'2N,2N-1V I 



J_ / 2V3(,,^(,) + ,^(,)) ^ \ 

7r2e2 \^ ^ iV2/3 ^ 7V2/3 ^ ^^'^ 'j ' 



/,(V2) / . ^1 ;r2e2 \ iV2/3 iV2/3 

"'2N-l,2N-2\' 'i-J \ / 

It follows that the integrands of ^2^^^ and l2^^ agree up to the order 0{N). According to (1.45) 
s = 2^1^ X + 0{N'~'^^^). Inserting the formula (2.18) into either integrals I^^^ or l2^\ we obtain 

+^ + 0{N-')^^drdu (2.19) 



/ / 25/3iv4/3r'(24/3x) + O(A-) dr dtx, 



where I2 is either of the integrals ^2^^^ and l2^\ If we write u = 1 — yN 2/3^ then we can write 
the integral (2.19) in terms of the variables x and y: 

h= / (25/3Ar4/3r'(24/3^) + 0{N))N-^'^dx dy 

^ ^ (2.20) 



= / / (25/3T'(2^/3a;) + ©(iV-^/^)) 



dx dy. 



If we let 1/3 < e < 2/3 then, after integrating, the error term goes to zero as A'^ — t- 00, provided that 
it is uniform for large x. This uniformity follows from the Riemann-Hilbert analysis, as discussed 
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in section 5. It follows that this integral has a limit as N ^ oo, which is 

/ / 2^/^Ti2^/^x)dxdy = - 2^/^T{2'^/^y)dy = T{x)dx. (2.21) 

JL Jy Jl 2 J24/3L 

Combining (2.16), (2.17), and (2.21), we thus find that, for a = 1 - LN''^/^, 

FN{a) = Io + h + l2 = -- T2N+i{x)dx + 0{N~^) 

T T^"- / X (2.22) 

= -- / R{x) + q{x) dx + 0{N-^) , 0<S< 1/3. 
2 724/3L V / 

A simple change of variables and exponentiation gives (1.25), the result of Theorem 1.1. The proof 
of Theorem 1.3 is very similar, and we present it in the next section. 

3. Proof of Theorem 1.3 
A rescaling of (1.38) using (1.31) gives the formula 

.2\ 2 



oa V2/ \hn-i,n-2[^-a) rin,n-i{a) 



(3.1) 



4"n-l(«) 

The small r behavior of the function Z^^^^\r) is described in the following lemma. 
Lemma 3.1. As r Q, 

4.0P.,,,,^Mr)_^,„,(_L),e-'.0(5). (3.2, 

where F^^^^ is defined in (1.56). 

The proof of this lemma is in Appendix A. Consider now Fn{a) for a close to 1. According to 
(3.1) and (3.2), we have 

Fr^^^)(a) = + - ilog (^) + r rG(^oPE)^^^^^^^ 

= f!^^^) - i log f ^"j + + 4^°^^) + + 

Zi \ ThTT J Z 

^' 2r2 4V /in-i,n-2(e-r)^+ K,n-l{r) 

(3.4) 



where 



and 



jiDOPE) ^ rG(^o^^\r)drdu, ^^^^ = T / G(^^^^\r) dr dn , 

Jo Jo Jl-n-^ Jo 

r r G('''^''^\r)drdu, 

Jl-n-^ Jl-n-^ 



j{DOPE) 
-'2 



(3.5) 



16 KARL LIECHTY 



for some 0<£<2/3. In the regime r < 1 — n for some < 5 < 2/3, the recurrence coefficients 
^rTfe are exponentially small in n, as shown in Appendix B. Therefore, inserting the asymptotics 
(2.14) into the integrals j^^^^^^ and j^^*^^^^^ tj^q £^(1 that both j^^*^^^) and /(^*^^^) are 0{n~^) 
as n — > oo. 

We now evaluate l2^^^^\ We scale o as a = 1 — Ln~^/^and write r = 1 — xn"^/^, so that as r 
varies from 1 — , x varies from r?l^~^ to L. Applying the asymptotics (1.50), we find that in 
this regime, 

i^A^:lM^) + = (^(^) - '/(^+)) + 0{n--l\ (3.6) 



which, using (1.53), is 0(n-2/3). From (1.48), we have that -^^P^ = 0(1), and thus 



^'^(e+4t,nK+) + <^i(r))' = 0{n-'l\ (3.7) 



Applying the asymptotics (1.55) the integral i^*^^^^ can thus be written as 

2i/3 



r r [-i7^(r;(.)+r;+i(.))+o(n-^) 



dr du. (3.8) 



According to (1.45), s = 2^/^x + 0{n ^^^), and we therefore have 

-l-Ln-2/3 



Writing u=l — yn and taking 1/2 < e < 2/3, this becomes 

2/3-e 2/3-e 



/•n2/3-e 

Jl 



It follows that 



= -^^og (^^2(2^/3^)^ +0(n-i- 



(3.9) 



(^!«)+0(n-V3-e)).. (3.10) 



(3.11) 



from which (1.57) follows immediately. 

The rest of the paper is dedicated to the proof of Proposition 1.2, which is based on steepest 
descent analysis of a discrete Riemann-Hilbert problem. 
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4. RiEMANN-HlLBERT ANALYSIS 

4.1. Equilibrium Measure and the 5- function. The equilibrium measure associated with the 
weiK ht e-"^(^) is the unique measure which minimizes the functional 

H{iy) = J j \og^-^^dv{x)dv(y) + j V{x)du{x), (4.1) 

over the set of probability measures on M. In the case that V{x) is given by (1.35), it is well known 
that the solution to this equilibrium problem is supported on the interval [— ttv^^ 
interval it has a density 

vra 

Y V TT^a 



dv^ix) = — \ a;2 dx. (4.2) 



Let us denote the density p{x). Clearly p{x) has its maximum value at .t = and p(0) = ^/a. The 
critical value of the parameter a, for which p(x) attains the upper constraint, is Oc = 1. 

We define the ^-function associated with these orthogonal polynomials as the log transform of 
the equilibrium measure: 

rb 



/ 2 
log{z — x)p{x)dx, b= — 1=, (4.3) 
-b TT-v/a 

where we take the principal branch for the logarithm. This function satisfies the following proper- 
ties: 

(1) g(z) is analytic in C \ (—00,6]. 

(2) For large z, 

00 „(, 2j 

g{z) = log z-J2^, 92j= ^dMx). (4.4) 

j=l ^ J-b J 

(3) The Euler-Lagrange variational conditions for the equilibrium problem (4.1) are 

r = v{x) + 1 for xe [-b, b] 
'^^'^^'-^'^\<Vix) + l for xeR\[-b,b], ^'-'^ 

where and g- refer to the limiting values from the upper and lower half-planes, respec- 
tively, and Z G M is the Lagrange multiplier. 

(4) The function 

G{x) = g+{x) - g^{x) (4.6) 
is pure imaginary for all real x, and 

G{x) = 2m ! p{s) ds. (4.7) 

J X 

(5) Prom (4.5) and (4.7) we obtain that 

2g±{x) = V{x) + l±2m j p{s)ds for x G [-6, b]. (4.8) 

(6) Also from (4.7), we get that G{x) is real analytic on the set (—6, b). We can therefore extend 
G into a complex neighborhood of (—6, 6), and the Cauchy-Riemann equations imply that 

dG{x + iy) 



dy 



= 2'Kp{x) > 0. (4.9) 

2/=0 
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Prom (4.5) we have that 

G{x) = 2g+{x)-V{x)-l = -[2g^{x)-V{x)-l], xe[-b,b]. (4.10) 
The value of the Lagrange multiplier is given by the equation 

1 



TT^ae 



(4.11) 



4.2. Interpolation Problem. The orthogonal polynomials (1.27) are encoded in the following 
interpolation problem (IP). Por a given n = 0, 1, . . ., find a 2 x 2 matrix- valued function Fn{z) = 
(P„(z)y )i<jj<2 with the following properties: 

(1) Analyticity: Fniz) is an analytic function of z for z £ C \ Ln^a- 

(2) Residues at poles: At each node x G L^^ai the elements P„(z)ii and Pn(^;)2i of the matrix 
P„(2;) are analytic functions of and the elements Pn(-z)i2 and Pn(^;)22 have a simple pole 
with the residues, 

Res Vn{z)j2 = Wn{x)Fn{x)jU 3 = 1, 2. (4.12) 

z=x 

(3) Asymptotics at infinity: There exists a function r{x) > on Ln,a such that 

lim r(x) = 0, (4.13) 

and such that as z — t- cx), Pn(z) admits the asymptotic expansion, 



p„(z)~(/+^ + ^ + ...j(; ^rj, ^Gc\ 



U D{^,r{^)) 

,OL 



(4.14) 



where D{x,r{x)) denotes a disk of radius r{x) > centered at x and / is the identity 
matrix. 

It is not difficult to see (see [3]) that the IP has a unique solution, which is 

where the Cauchy transformation C is defined by the formula, 

^(/)(-)= E (4-16) 

,OL 

Because of the orthogonality condition, as z — )• oo, 

-C{WnPn){z)= > > ^n{x) y -r^i = + > (4.17) 

n ^-^ niz — x) ^-^ n ^-^ z^+'- z^+'- ^-^ z^ 

xeLn,c, xeLn,a J=0 J=n+2 

for some constants aj, which justifies asymptotic expansion (4.14), and we have that 

hl^l = [Pi]i2, (/iS-i)"' = [Pi]2i. (4.18) 
Furthermore, the recurrence coefficient A^^_^ is given by (see [3], [7]) 

[P1J21 

and the Christoffel-Darboux kernel (1.41) is given by 

-nV{x)/2 -nV{y)/2 /-.X 
^"^"^'^^^ ^(^^^) (° 1) Pn(ar)-^Pn(y) y . (4.20) 
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The asymptotic analysis of this IP follows the steepest descent method of Dcift-Zhou, the plan 
of which is as follows. Wc first convert the IP to a Riemann-Hilbert problem (RHP), where the 
condition on poles and residues is replaced by a jump condition on some contours in C. Then 
we perform a series of explicit transformations which convert the RHP to one with jumps which 
approach the identity matrix as n oo. This small norm problem can be solved by a series 
of perturbation theory. We can then recover the orthogonal polynomials encoded in the IP by 
inverting the explicit transformations which led us to the small norm problem. 

4.3. Reduction of IP to RHP. We now reduce the interpolation problem to a Riemann-Hilbert 
problem. Introduce the function 

sin(n7r2; + air) 



mr 



Notice that 

n(xjk) = 0, Il'{xk) = exp (imrxk + ina) = (—1)'^, for Xk = 
Introduce the upper triangular matrices, 



k — a 



n 



' nn{z) ^ 



(4.21) 
(4.22) 

(4.23) 



and the lower triangular matrices, 

U{z)-^ 
± ~ \ n ^ ^±i7r(nz+a) n(z) 



D 



Define the matrix-valued functions, 









Uiz) 



U{z)Wn(z) 



^±iTr{nz+a) 



and 



Rj, = Pniz) X 



D"(z) when Imz > 
when Imz < 0, 



D'_(z), 



when Imz > 
when Imz < 0. 



(4.24) 



(4.25) 



(4.26) 



The functions Il^{z), tin{z) are meromorphic on the closed upper and lower complex planes and 
they are two-valued on the real axis. Their possible poles are located on the lattice -Ln,a- A-S shown 
in [7], in fact they do not have any poles at all. 

Consider the regions Q,^ and fij shown in Figure 2, and the contour S which bounds these 
regions. These regions lie entirely within an e-strip of the real line, and the regions are small 
sectors above and below the origin, respectively, within this strip. We make the transformation 

KRl{z)K-\ for zenj, 
Kniz) = { KR^(z)K-i, for z e (4.27) 
KP„(z)K~^, otherwise. 



where K 



'1 
^0 -2i7r^ 

Let us denote by 7± the part of the contour S which is the boundary between the region and 
the region rzj. Then the region $7^ is bounded by the contour 7-1- and a small segment on which 
Imz = =be. Let us denote these segments a±. 
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Figure 2. The contour S dividing an e- neighborhood of the real hne into the 
regions fl^ and fij. 



4.4. First transformation of the RHP. Define the matrix function Tn{z) as follows from the 
equation 

R„(^) = e'5^'^3T„(z)e'^(^(")-5)-3^ (4.28) 

where I is the Lagrange multiplier, the function g{z) is described in section 4.1, and o"3 = 

is the third Pauli matrix. Then Tn{z) satisfies the following Riemann-Hilbert Problem: 

(1) T„(z) is analytic in C \ S. 

(2) Tn+iz) = Tn-iz)jT{z) for z G S, where 



jriz) = { 



(3) As z — )• oo. 



4.5. Second transformation of the RHP. Introduce the matrices 



(4.29) 



(4.30) 



A+{z) 



_J_^—in{nz+a) a 

2mTi^ . ^ 1 \ (7) 

-2n7rie^"("^+") / ' ' 



1 ATrlnz+a) 
2mri 









2n7rie-^'^("^+")^ 

(4.31) 

We now make the transformation 

(Tn{z)U{z)-^ for zG{{-b,b)x{0,ie)}nnX 
Tn{z)j-{z) for ze{(-6,6) X (0,-ie)}nJ]Y 
T„(z)A±(z) for 2 G 
^ T„(z) otherwise. 

This function satisfies a RHP similar to T, but jumps now occur on a new contour S5, which is 
obtained from S by adding the segments (—6 — ie, —b + ie) and (6 — ie,b + ie), see Figure 3. On 
this contour, the function S„(z) satisfies the jump condition 

S+iz) = S^{z)js{z), (4.33) 
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Figure 3. The contour Sg. 



where 



js{z) 



1 
-1 

1 ^n{g+{z)+g-iz)-l-V{z)) 

1 



when z G (-6, b) \ {0} 

when z E M \ (-6, 6) 

when z G {(—6, 6) it ie} \ a± 



^1 _ g±2i7r(n2+Q)-j-l _|__ 



p±nG(z) 



_ g±27ri(n2-|-a)^-l 







TnG(2)) 



(4.34) 



1 _ ^±2ni{nz+a) ) ^hen z G a± 



]^ ^^±nG{z)±2TTi(nz+a) 



when -2 G 7± 



i_eT2.i(n^+-) z G UM\ (-6,6)} ±ie I . 





1 ±' 


By formula (4.5) for the G-function and the upper constraint on the density p we obtain that, 
for sufficiently small e > and x G (—6, —e) U (e, 6), 

< ±ReG{x ± ie) = 2Trep{x) + 0(6^) < 27re + 0(6^). (4.35) 

Combined with property (4.7) of the ^-function, this implies that all jumps on horizontal segments 
are exponentially close to the identity matrix, provided that they are bounded away from the 
segment (—6, 6). 

4.6. Model RHP. The model RHP appears when we drop in the jump matrix jsiz) the terms 
that vanish as n — ?• oo: 

(1) M{z) is analytic in C \ [-6, 6]. 

(2) M+{z) = M_{z)jM{z) for z G [-6,6], where 



jMiz) 



1' 

-1 0/ • 



(3) As z — )• oo. 



, , , Ml M2 

z 



This problem has the unique solution (see e.g., [16]) 

M{z) -- 

where 

l{z 



j{z)+j-Hz) 7(2)-7-l(2) ' 

2 -2i 

-l{z)--/-\z) 7(2)+7-l(z) 

2i 2 / 



z + b 



z-b 



1/4 



(4.36) 
(4.37) 

(4.38) 

(4.39) 
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with a cut on [—b, b], taking the branch such that 7(2;) ~ 1 as z — >■ 00. 

4.7. Parametrix at band-void edge points. We now consider small disks D{b,s) and D(— 6, e) 
around the endpoints of the support of the equilibrium measure. Denote 

D = D{b,£)iJD{-b,£). (4.40) 

We seek a local parametrix U(z) defined on D such that 

(1) 

\]{z) is analytic ouD\ S5. (4.41) 

(2) 

\5+{z) = \5-{z)js{z) for zGDnS^. (4.42) 

(3) 

U(z) = m{z) {I + O(n-i)) uniformly for z G dD. (4.43) 
The solution to the problem is given in [16] (see also [7]), and we do not repeat it here. 

4.8. The Riemann-Hilbert problem associated with the Painleve II equation. In section 

4.9. we will discuss the local analysis to our Riemann-Hilbert problem near the origin, which is the 
point at which the equilibrium measure attains the upper constraint. The solution will be given in 
terms of a well known problem in integrable systems, which we discuss now. For a more complete 
description of this problem, see the book [23]. 

Let *(C) be the 2x2 matrix- valued solution to the differential equation 

dC^^^^I \ 4Cq-2ir 4iC^ + i{s + 2q^) J ^ ' ' ^^'^^^ 

It is known that there exist solutions *j {Q defined in each of the six sectors 

5, = |cGC:^:^7r<argC<^^7r|, (4.45) 

such that for j G {1, 2, 3, 4, 5, 6}, as ^ — )■ 00, 

*^-(z)e^(tC'+<)^3 ^ J + 0(C~^), (4.46) 
and on the ray Tj = : arg ^ = ^^^tt} we have the jump condition 

*,+i(C) = *,-(C)^,-, (4.47) 

where 

^^={a, !) ^^ = (0 ?) 

The numbers Uj are called the Stokes multipliers and satisfy the relations Oj+s = aj, and 010203 + 

Ol + «2 + «3 = 0. 

If the parameters s, q, and r are chosen such that q, as a function of s, is the Hastings-McLeod 
solution to the Painleve II equation and r = q'{s), then the Stokes multipliers become 

ai = l, 02 = 0, 03 = -1. (4.49) 

It is known that the Hastings-McLeod solution to the Painleve II equation is a meromorphic function 
whose (infinitely many) poles are all located away from the real line. We will consider only real s, 
and thus we can always choose the parameters r, q, and s in such a way. Because q{s) behaves like 
the Airy function for large s > 0, in fact the error in (4.46) is 

O I — - — I , as C ^ c)o, +00. (4.50) 
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Figure 4. The contour T, which is composed of the rays 7j, j = 1, 2, 3, 4. 



Consider the oriented contour F in the C-plane depicted in Figure 4, which is made up of the 
four rays 71,72,73, and 74, where 



71 



73 



{CGC:argC=^} 



C G C : argC = — 





72 



74 



C G C : arg C = — 

D 



llvr ' 

CGC:argC = — !> • 



(4.51) 



The function ^'(C) which solves differential equation (4.44) and corresponds to the Hastings- 
McLeod solution to the Painleve II equation satisfies the following Riemann-Hilbert problem. 

(1) *(C) is analytic for C G C \ F. 

(2) For C G F, ^{C) satisfies the jump condition 



*+(C) = *-(C)j*(C), 



where 



MO 



1 
-1 1 

1 -1 
1 



for C G 71 U 72 



for C G 73 U 74 , 



(4.52) 



(4.53) 



(3) As C ^ 00, 



*(() = (/ + 0(C-^))e-^(t^'+^^)"3 



Introduce the parameter ^ G M, and let *2(C; 9) be defined via 'J'(C) as 



*2(C;s,0) 



Vi ImC>0 



(4.54) 



(4.55) 



where ai = and (T3 = are Pauli matrices. Then the function ^2{C]s,9) satisfies 

the Riemann-Hilbert problem 

(1) *2(C; 0) is analytic for C G C \ (F U M) . 

(2) For C G F, *2(C; S;^) satisfies the jump condition 



*2+(C;5,^) = *2-(C;5,^)i2(C;s,^), 



(4.56) 
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where 



1 _e2i(|C^+<-e)cT3 

1 

1 g-2i(|c^+<-e)<73 
1 



C e 71 u 72 

C G 73 U 74 



(4.57) 



1 
-1 



(3) As C ^ oo, 



*2(C;s,^) 



^ {I + 0{C')) J) ImC>0 

(/ + 0(C-^)) ImC<0. 
Finally we note, as in [9], the formula for the entries of the matrix 



<fi(C;s) 



C e Si u S4, 



(4.58) 



(4.59) 



where $1,2 (C; s) are defined in (1.61). 

4.9. Parametrix at the origin. We seek a local parametrix U(z) defined on D{0,e) such that 

(1) U(z) is analytic for C G D{0,e) \ S^. 

(2) For C G r, *2(C) satisfies the jump condition 

V+{z) = V.{z)js{z) for zGi5(0,e)nSs. (4.60) 

(3) As n — >■ 00, 

U(z) = M(z)(/ + 0(n-i/3)) uniformly for z G 9i:»(0,e). (4.61) 
Let us recall the jumps js on the contours 7+, 7_, and M close to the origin: 

'j^ ^nG{z)+2ni(nz+a) 



1 

2 g— nG(z)— 27ri(nz+a) 

1 

1 



z e 7+ 

2; G 7- 



(4.62) 



-1 



z G 



These jumps are similar to the jumps given in (4.57). We thus seek C,{z;t), s{t), and 9{t) which 
solve the equation 



2i ( ^C{z; af + s{a)C{z; a) - e{a) ) = n{G{z) + 2Triz) + 27ria, 



or equivalently 



1 ^, ,0 s(a) ^, , Qia] 
-C(^;a)3 + ^C(^;t)-^ 



-(G(z) + 27rzz) + — 



(4.63) 



(4.64) 



As shown in [8], there is a unique solution to this equation which is regular at the origin, and it is 
defined in terms of the stationary points of the right side of (4.64). If we denote 



f{z;a) = -iG{z) + 27nz) + ^ = ^ 



1 + 2^-2 



/ p{x)da 
Jo 



+ 



4 ' 



(4.65) 
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then the zeroes of the function f'{z;a) (here ' means differentiation with respect to z) are at the 
stationary points 

2 2 

zi{a) = — va — 1, Z2{a) = va — 1. (4.66) 

(4.67) 
(4.68) 
(4.69) 
(4.70) 



ira vra 
Notice that 2;i(l) = 2:2(1) = 0. We then have that (see [8]) 

e{a) = -2(/(zi; a) + f{z2; a)) , s{af = -36{f{z2;a) - /(zi; a))'. 
Notice that, as p{x) is an even function of x, we have 

mrz 



f{z; a) + f(-z- a) = ^ + f{z; a) - f{-z; a) 



It follows that 



e{a) 



riTT 



— ira, s{a)'^ = —dir^n^ ( zi 



j-zi 

— I p{x)dx 



The possible zeroes of C,{z; a) are the solutions to the equation 
or equivalently, 



, rivr 7ra 
f{z-a) = — + —, 



= I p{x)dx. 
Jo 



(4.71) 



The only solution to this equation, and thus the only possible zero of ({z), is at 2; = 0. Indeed, it 
is not difficult to see that ^(0; a) = 0, and we thus have 



We now take 



s(a) 



tS{z;s,0) = < 



s{a) 



s{a) 



(4.72) 



J -.v,v-.,-,-. ^4^3^ 

{M{z)^2iCiz);s,e) lmz<0. 
This function is analytic for z G D{0,£) \ S5. For z G S5, it satisfies the jump conditions 



where 



ju{z) 



lJ+{z;s,e) = V^{z;s,e)ju{z) 



1 

J g-2i(|C(2)^+<W-e)a3 

1 

1 
-1 



(4.74) 



z G 7+ 



(4.75) 



C G M. 

Let us check the large n behavior of the function \i{z). As n ^ 00, we have 

miz 



\J{z;s,e) = { 



\ -^){I + 0{n-^/'))[^^ J) lm.>0 
M.{z){I + 0{n-^/^)) Im2<0. 

m{z){I + o{n-^'^)) . 



(4.76) 



It follows that we may take our local solution to be U(2;) = \J{z\s,9), where s = s{a) and 9 are 
given in (4.69). 
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X„(z) 



(4.77) 



4.10. The third and final transformation of the RHP. We now consider the contour Ex, 
which consists of the circles dD{—b,e), dD{b,e), and dD{0,e), all oriented counterclockwise, to- 
gether with the parts of S5 \ [—6,6] which lie outside of the disks D{—b,e), D{b,e), and D(0, e). 
Let 

' Sn{z)M{z)~'^ for z outside the disks D{-b,e), D{b,£), D{0,e) 
Sniz)\J{z)~^ for z inside the disks £>(-&,£), D{b,e), D{0,e). 
Then X„(^;) solves the following RHP: 

(1) 'X.n{z) is analytic on C \ Sx- 

(2) X.n{z) has the jump properties 



where 



(3) As z ^ 00, 



jxiz) 



M(z)U(z)-^ for z on the circles 
M{z)jsM{z)~^ otherwise. 



X„(.)^/ + ^ + ^ + ... 



Additionally, we have that jx{z) is uniformly close to the identity in the following sense: 

I + 0{n~^) uniformly on the circles dD{—b,e), dD{b,e) 
jx{z) = \ I + 0{n~^^^) uniformly on the circle dD{0,s) 
I + O (e-'^(^)" ) on the rest of Ex , 



(4.78) 
(4.79) 

(4.80) 
(4.81) 



where C(z) is a positive function bounded away from zero and with sufficient growth at infinity so 
that e-^(^) G L^{^x)- 



If we set 
then (4.81) becomes 



fxiz)=Jx{z)-I, 

0{n^^) uniformly on the circles (9 Z)(— 6, e), dD{b,e) 
0(n^^/'^) uniformly on the circle dD{0,s) 
^ 0(e-C(^)n) on the rest of Ex- 



(4.82) 



(4.83) 



The solution to this small norm problem is given by a series of perturbation theory. Namely, 
define the functions Vk recursively as 



Vk{z) = - 



2m 



Vk-i{u)jx{u) 



du, vo{z) = I, 



(4.84) 



where Z- means that the integration takes place on the minus-side of the contour. The solution is 
then 



k=l 



where 



^n{z) = I + ^ ^n,k{z) , 

Vk-i{u)j%{u) 



X„ u(z) = 



du. 



z — u 



(4.85) 
(4.86) 
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In particular, this implies that 
uniformly for z G C \ T,x ■ 



5. Evaluation of Xi in the critical case 
The function *(C) satisfies, as C ^ oo. 



e-i{|C^+<)^3^(^)-i = ^ 



where 



A = -R{s)a3 + q{s)aia3 , B = {q{sf - R{sf)l - 2{q'{s) + q{s)R{s))ai , 
q{s) is the Hastings-McLeod solution to Painleve II, and 



/•oo 

R{s)= / 

J s 



see [19]. For u G dD{0,e), we have that 



jxiu) = 



-1 



Imu > 



Prom (5.1) and (5.4), we see that 

M{u)aie' 



^ A + TTT^B + O(n-i) ) e-^^^ViM(«)-i 



Im« > 



2i<(u) 8((u)2 



Im u < . 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



Notice that, according to (4.50), the error term in (5.5) is uniform for s G [so,oo) for any sq G M. 
As can be seen in Appendix 2.3, the jump matrices on the contours D{±b, e) are uniformly close to 
the identity for a close to 1. Thus if wc write a = 1 — xn~'^^^, using (1.45) and (4.85), any error we 
compute in the large n expansion of X„ will be uniform for x G [L, oo) for any L G M. Multiplying 
out the above expression, we get 



(5.6) 
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R{s) 



7(i()^ + 7(n) ^ ]cT3 + ^^(j{uy^ -jiu)''^ lo-i + (-l)"cos(27ra)g(s)o-3£Ji 



i-iy 



2i 



sin(27ra) ( j{u)'^ — j{u) ^ I q{s)a3 



i-rr 

2i 



sin(27rQ;) j^u) + j{u) q{s)ai , 



R2{u) = {q{sf - R{sf)l - cos(27ra) (^j{uf - 7 {q'is) + q{s)R{s))a3 

- i-ir cos(27ra) (^7 W + tW"') {q'{s) + qis)R{s))ai 

- 2z(-l)"sin(27ra)(g'(s) + g(s)ii;(s)) (73(71. 

(5.7) 

anci 7('u)^ is the analytic continuation of j{u)'^ from the upper half plane. The Taylor expansion 
of 7('u)^ about the origin is 



7(n)2 

which in turn gives that 



j{uf - 7^-' = -i[2 + - + 0{z')j , ^{uf + 7(«)-' = 
Expanding each term in (5.6) about the origin gives 

.0 , , 1 (Rm ^ tr,. ^ (R'm r(o)^i(o) 



Let us write this expansion as 
1 fA 



n 



1/3 



In particular, 



, n^/^Ri{0) . n^/^R[{Q) ^ n^/^R'^i^) ^ n 



2<'(0) ' 2iC'(0) 
Let us evaluate X.n,i{z). We have 



8C'(0)2 



1 



1(2;) = 



U 



du + 0{n~^). 



ldD{0,e) 

This integral can be evaluated via a residue calculation. Indeed, for z G C \ D(0,e) 



1 



j\{u)du ^ _1 
z 



'dD{0,e) 

In particular, this gives that, for z G 

1 

2; n 



+ 0(n-i). 



^1(2) 



1/3 „2/3 



+ 0{n-'), 



(5.8) 



(5.9) 



(5.10) 



^_ ^ ^ _ r(0)i?2(0) , \ , 

^8C'(0)2V n2 + + 2 3C'(0) ^ 

expansion as 

(^ + Ao + A^u + ---)+47^(^ + ^ + Bo + ---) +0{n-'). (5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 
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y^n,2{z) = Tr- 
im Jq 



dD{0,e) 



and thus 

u — z 

1 

27ri JaDio,e) (u - z)u 
For z G C \ -D(0, e) this is evaluated as 



n 



1/3 



(5.16) 



n 



2/3 



fx{u)du + O(n-i). 



+ 0(n-i) 



(5.17) 



This can also be taken as a formula for V2{z) for z G S^. 

Let [X„^jfc]j be the coefficient of the term in the expansion of X.n^k{z) at z = oo, so that 



Xn,fc(^) = 

k. 

From (5.15) and (5.17), we see that 

Adding these together gives that 

Xi = [X„,i]i + [X„,2]i + 0(n-^) 



[Xn,fc]l _|_ [Xn,fc]2 _|_ 



fx. 



n,2jl 



^-1^0 

n2/3 



+ 0(n-^) 



+ 



1 



j^l/3 77,2/3 
1 



ii?i(0) - 



-B^i + A_iAo) + 0{n~^ 
R'^iO) Ri{0)R[{0) 



4C'(0) 



2C'(0) 

Also from (5.15) and (5.17), we have 

X2 = [X„,i]2 + [X„,2]2 + 0{n~^ 

B-2 A^ 1 1 
j^2/3 77,2/3 ^ ^ 



2C'(0) 



0{n-^). 



1 



8C'(0)2 



[i22(0) + 2i?i(0)2] +0(n-^). 



Notice that 



^i(O) 

it;i(o) 



-i?(s)o-i + (-1)" cos(27ra)g(s)cr3cri + i(-l)" sin(27ra)g(s)c73, 
-«^3 + ^ -TT^ sin(27ra)g(s)(7i, 



2i " 2 
i?2(0) = (g(s)2-i?(s)2)/-2(-l)"cos(27ra)(g'(s) + g(s)i?(s))(73 

- 2z(-l)" sin(27ra)(g'(s) + q{s)R{s))azai , 

R'^{(}) = {-l)''^^^^/^cos{2^^a){q'{s) + q{s)R{s))ai. 



(5.18) 
(5.19) 



(5.20) 



(5.21) 



(5.22) 



6. Proof of Proposition 1.2 and Theorem 1.4 

The quantities /i„ „ and /in,n-i are encoded in the matrix Pi described in (4.14). According to 
(4.27), (4.28), and (4.77), 

p -K-i/" [Xi]ii + [Mi]n ([Xi]i2 + [Mi]i2)e"'\ . . 

^^-"^ V([Xl]21 + [Mi]2l)e-«' [Xi]22 + [Mi]22 J ^ ^ 
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It follows from (4.18) that 



hn,n = -27ri([Xi]i2 + [Mi]i2)e 



nl 



-2m 



and 



^/a \7r^oe 
1 



(l - [X.l\l2^T^^/a) , 



(6.2) 



n,n—l 



-2m 
1 

~2m 



([Xi]2i + [Mi]2i)e- 



/ 2 

[TT ae) 



[Xi]21 



-nl 

lb 
2" 



(6.3) 



^ ^ n (vr^ae)" (l + [Xi]2i7riv^) . 

2^Ja'K^ 



According to (5.20) and (5.22), 

22/3 



[Xi]i2 m\fa 



and 



[Xi]2i m\fa ■ 



nV3 

+ 

22/3 



R- (-l)"cos(27ra)g 
2i/3 

2/3 



n 



(-I)'* cos(27ra)(g' + 2qR) - + sm^{2'Ka)^ + 0(n"^), 



(6.4) 



n 



1/3 



i? + (-l)"cos(27ra)g 

2V3 / \ 
+ ( (-1)" cos(27ra)(Q' + 2qR) + R^ - q^ sin2(27ra) J + 0{n 



(6.5) 



where the functions q and R written with no argument refer to those functions evaluated at s. This 
proves equation (1.48). 

Also from (4.27), (4.28), and (4.77) we get 



[P2]21 . [X2]21 + [M2]21 + [XiMi]21 

— = FiFl , nv,r 1 L-^lJll ~ [-I-Vlljll. 



[Pl]21 ""^^ [Xi]21 + [Mi]21 ^"^^^^ ^^-^^ 

Asymptotic evaluation of this expression from (5.20) and (5.21), in light of (4.19), proves (1.50). 
Thus Proposition 1.2 is proved. 

We now turn to the proof of Theorem 1.4. Consider z G D(0, e) H JlJ. In this region we have 

'K-ie'^"3X„(z)M(z)aie*^"^^*(C(2))e*(^(")'+^^(^'*)-'^)"Vij+(z)e"(5(")-'/2)-3 
xK(D"(2;))"^ lmz>0 

Pn(^) ={ ^ +^ " ^gj) 

K-ie'^'^^^X„(z)M(z)e^^'^V3*(C(^);s)e^(^(^) +^^(^'*)-^)'^V3j_(z)-^e'*(^^(^)-^/2) 

xK(D!l(z))-^ lmz<0. 

Of course C,{z) and s depend on the parameter a, but we suppress the notation here. Make the 
scaling 

c = 7r2-V3, a = l-Ln-^/K (6.8) 



Notice that in this scaling 



y 



an 



1/3' 



C(x) = u + 0(n-2/3), C(y) = ^^ + 0(n-2/3), s = 22/3l + 0(n-2/3)^ 



(6.9) 
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and that 



It then fohows that 



s(l) = Soo= Hm sit) = 



(6.10) 



X I + O 



u — V 



n 



1/3 



X e"(s(2/)-'/2V3KD't(y)"^ lmx,y>0 
D!i(x)K-ie-"(f(^')-'/2)a3^-_^^)^gg-K|CW«+<(a:)-0)^3^^^^^).g)-i 

\ ~ _i(4^(y)3+sC(y)-e)cT3^ /.A-1 



(6.11) 



X [I + O 



n 



1/3 



*(C(y);s)e^(3< 



'<^3j-(y)" 



Taking hmits from either the upper or lower half planes, we find that for x and y real and in -D(0, e), 



Kn{x,y) 



-nV{x)/2„-nV{y)/2 



n{x - y) 

gj^l/3 

2imr{u — v 

gj^l/3 



nG(x) 

-e 2 



n — v 



n 



1/3 



(0 l)P4x)-iPn(y)(^Q^ 

e-=^) e-^(t^(-)'+<(-)-^)-3*(C(x); s)-i 
*(C(y);s)e^(l< 



/ nG(i<) ■ 
1 "G(h) 



(6.12) 



e 2 



2imr{u — v) 



g-m(nj/+a)J 



In order for this kernel to make sense, we must choose x and y to lie on the lattice If we take 

u ^ V and 



X = 



kn- a 



n 



y = ; kn,mne^, 

n 



(6.13) 



u — V 



n 



1/3 



then the above expression becomes 

/'_-|^fc„+m„ 1/3 / / 

According to (4.59), 

*^^'^)0) = (I^(c!2))' i)*-HC;^) = (-MC;^) $i(C;^)), 

thus (6.14) imphes (1.65). Notice that the presence of the factors e='=*'^("'^+") in (6.12) plays little 
role in the limiting values of the off diagonal terms. To compute the diagonal terms, we must take 
a limit of (6.12) as u ^ v, and these factors do indeed play a role. In taking this limit, we must 
take into account the fact that defi^ = 1, and we obtain (1.67), which proves Theorem 1.4. 



*(C(y);s)L . (6.14) 



(6.15) 
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Appendix A. Proof of Lemmas 2.1 and 3.1 

The proofs of Lemmas 2.1 and 3.1 are based on the following estimate of the rate of convergence 
of a Riemann sum, which is slightly sharper than the a priori rate of 0{e). 

Lemma A.l. Let f{xi, . . . ,Xn) be an analytic function of n variables. Let the functions A^ixi, ... ,Xn) 
be defined recursively via 



d 



Ak{xi,...,Xn) = ^-Q^Ak-l{xi, . . . ,Xn) , Aq = f. 



(A.l) 



Suppose that 

■■■ Ai{xi,...,Xn)dxi---dXn= ■■ A2{xi, . . . , Xn)dxi ■ ■ ■ dXn = 0. 

Then as e ^ 0, 

/•••/ f{xi,...,Xn)dxi---dxn-e'' V = 0(£^). 



(A.2) 



(A.3) 



Proof. A multi-integral my be estimated by writing 



/r rxi+e rxn+e 

■■■ f{xi,...,Xn)dXi---dXn= /J / f {ti, . . . ,tn)dti ■ ■ ■ dtn. (A.4) 

Expanding each integrand on the RHS as a Taylor series and integrating term by term, this becomes 



Xl,...,Xn&eZ 



where 

/ = fi,-^li ■ ■ ■ ) -^n) ) -^k = -^jixii ■ ■ ■ 1 Xji) ■ 

Thus the error in the Riemann sum is given by 

/•••/ fdx^...dxr.- Y: ^"/ = ^" E fT^ + T^' + #^' + 0(^') 



(A.5) 



(A.6) 



3^1 ,...,a;7xt 



(A.7) 



X\^...^Xji^£'^ 



X\ ,...,Xn 



If (A.2) holds, then by the same argument, 



^ ^i(xi, . . . ,x„) 
£" ^ A2(a;i,...,x„) 

a;i,...,a;ne£Z 



3^1 7 • • -i^^nSS^ 



2 



X\,...^Xji^£'^ 



(A.8) 



a;i,...a:nS£Z 
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and (A. 7) can be written as 



/•••/ fixi,...,Xn)dxi---dXn-e"' y2 f{xi,...,Xn) 
•1 •''M" „ 



x\,...x„&e1 

(A.9) 



□ 



, 12 I ^ / 24 



Let us first apply this result to the proof of Lemma 3.1. By rescaling we can write (1.34) as 

xi,...,Xn&e{L-a} l<j<k<n 



n 



where 



We have here an exphcit prefactor times a Riemann sum for the function 

(A.12) 



n 

X] 



f(xi,...,Xn) = Yi i^k - Xjf exp l-^ 

^<j<k<n j=l 

which is exactly the integrand in (1.56). One may easily check then that 

Ai{xi,...,Xn) = -2{xi H ^Xn)f{xi,...,Xn), 

A2{xi, ...,Xn) = 2f{xi, . . .,Xn){2{xi H \- Xnf - Tl 



= -2 Ai{xi,. . . ,Xn){xi H \-Xn) + uf (xi, . . . , Xn) 



(A.13) 



It is easy to see that Ai has the symmetry Ai{xi,. . . ,Xn) = —Ai{—xi, . . . , —Xn), from which it 
follows that 

/ • • • / Ai{xi, Xn)dxi ■ ■ ■ dxn = 0. (A.14) 
It is a simple exercise to integrate by parts to see that 

■■■ A2{xi, Xn)dXi ■ ■ ■ dxn = 0. (A.15) 

J jRn 

It then follows that as r — > 0, 



(A.16) 



Taking logarithms gives (3.2). 
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We now prove Lemma 2.1 in the absorbing case. The proof in the reflecting case is nearly 
identical. Using symmetry about the origin and a rescaling of (1.15), we get 



maxftS^^Vi) <m] = 
o<t<i ^ ' J 

2N(N+1) 



N 



where 



( N )\ (A.17) 



e = (A.18) 
My/2 

We again have an explicit prefactor times a Riemann sum for the integral 

foo i-oo f ^ \ f ^ 1 

■ ■ ■ (A(x2)) ' ( n x] 1 exp I - ^ x) Uxi • • • dxN. (A.19) 

This integral is the partition function for the Laguerre unitary ensemble. Its value is known (see 
e.g., [24]), and it exactly cancels the prefactor, so that 

lim P| max bN{t) < M] =1. (A.20) 

M^oo \0<i<l J 

Lemma A.l also holds for multi- integrals over M" , that is if we replace M"" with and Z with N, 
and we can thus use Lemma A.l with 







N 








exp < 




> . 


(A.21) 













It is not difficult to see that in this case the condition (A. 2) is satisfied. Indeed, notice that for any 
j = l,2,...,iV, 

/ (^•^(^i'---'^-'^^) ^ -f{{xi, . . . ,Xj-i,0,Xj+i, . . . ,xn) = 0. (A.22) 



Furthermore, notice that 
d 



dxj 



/(2;i,...,Xiv) =exp|-^x|| ^ JJ P(a;i, . . . , x^v) (A.23) 

for some polynomial P. It follows that, for any j = 1,2, . . . , N, 

Ai{xi, . . .,Xj-i,0,Xj, . . .,xn) = 0, (A.24) 

and thus 

/ ■ ■ ■'^n)^ dxj = -Ai{xi,. . . ,Xj-i,Q,Xj,. . . ,xn) = (A.25) 

(A.22) and (A.25) imply (A. 2), and thus Lemma A.l applies. It follows that 

m^^^ i^Sv ^^(i) < = 1 + 0{e^) = 1 + 0(M-^). (A.26) 
In the scaling M = ^ this becomes, as a ^ 0, 

P ( max (t) < m) = 1 + 0{e') = 1 + O . (A.27) 
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Taking the logarithm proves Lemma 2.1. 



Appendix B. Proof of Lemma 2.3 
If the parameter a is such that 

-X 2 
a<l-n °, <d < -, 



(B.l) 



then by (1.45) and (4.50) the jump matrix for X„(2;) about the origin is exponentially small in n, 
and therefore the asymptotic expansion for X„(2;) comes from the jumps on the circles dD(zizb,£). 
We need to calculate this expansion up to an error of the order n~^. Instead of doing this directly, 
which is rather tedious, let us proceed by comparing our discrete orthogonal polynomials with their 
continuous brethren, the monic scaled Hermite polynomials {Pj'^\x)}'^q, for which we have exact 
formulas. These polynomials satisfy the orthogonality condition 



p\''\x)Pi'\x)e 



(B.2) 



The superscript (c) stands for continuous. The continuous orthogonal polynomials Pn^ can be char- 
acterized in terms of the following Riemann-Hilbert problem. We seek a matrix P^'^^ (z) satisfying 
the following properties. 

(1) V''n\z) is analytic on C \ M. 

(2) For any real x, 

1 wix) 



P(f)(x) = p('=)(:c)j,(x), jc(x)=(j ^(^^ 



(3) As z ^ GO, 



where Pj. , A; = 1, 2, ... , are some constant 2x2 matrices. 
This problem has the unique solution 

p(c)/ \ 1 r Pjf' (u) w„(u) du 



>(c) 



+ ... 



(B.3) 



(B.4) 



p(c)c 



X] 



u—z 



(B.5) 



The normalizing constants /in ■* can be found as 



-2m 



12 



p(c)" 






21 


2m 



(B.6) 



(c) . (c) 

We can make a series of transformations to P„ to arrive at a small norm problem. Define T„ ' 



from the equation 



and Sn ■* as 



T(f)(z)i+(z)-i for ze{{-b,b)x{0,is)} 
Si:\z) = {Tif>{z)j^{z) for ze{{-b,b)x{0,-te)} 
tI^\z) otherwise. 



(B.7) 



(B.8) 



36 



KARL LIECHTY 



(B.9) 



(c) 

Then the matrix ' can be defined as 

J Sl^\z)M{z)-^ for z outside the disks -D(-6,e), D{b,e) 

\ Sl^\z)\J{z)-^ for z inside the disks D{-b,e), D{b,e) . 

The jump matrices for 'Kn\z) are exponentially close to those for X„(z), and therefore, by (4.85) 
and (4.86), X.'n\z) and X niz) are exponentially close to eachother. We are interested in the off 
diagonal terms of the matrix X^*^^ , where 



X^\z) = I + ^ + 0{z-'). 



(B.IO) 



One can easily see that 



X 



12 



-nl 



12 



[M 



1 12 



X 



21 



>(c) 



:^(^2^e)"-^, 
2m Try a 



21 

27rz 



- [Mi; 



21 



(B.ll) 



+ 



The constants and are known exactly: 

(,) _ _ n\V2^ 
" "(v 

It follows that 



(V^7r)2"+i ' 



"n-l 



/i2.i(7r2ae)'^ ^V^ 
(n - l)!V27f 



XWI = ' (l (^Y ^' ^ 

- ^ Ji2 ^^^/a\ \nJ -v/27m/ 



Applying Stirlings formula, we find that 

(J- 



X 



21 



TT-y/a \ \eJ ^/n{n—l)\ 



(B.12) 



(B.13) 



xi'^)" 






12 


TT-s/a \ 


xS'^)" 




' ( 


21 


7^^/a \ 



+ 



139 



139 



21 vrTa V 12n 288n2 51840^3 



+ 0{n-^) 
+ 0(n-4) 



(B.14) 



Let us now return to the discrete system of orthogonal polynomials. The normalizing constants 
are given as 



(. 



1 



, I (l - [Xi]i27riv^) , 

^"•"-^ " (TT^ae)" (1 + [Xi]2i7r?V^) . 



(B.15) 



Since X„ and X^'^^ are exponentially close to eachother, we may use the above expansion for Xi, 
obtaining 



1 



1 1 

1 + T7^ + 



139 



12n 288n2 51840n3 



+ 0(n" 



2^7r2 



W'aeY 1 



1 1 139 ^, _4, 

I I u Oin 

12n 288n2 51840n'^ ^ 



(B.16) 



Let us also note that in this asymptotic regime, by a similar argument, wc find that the recurrence 
coefficients A^j^(a) are exponentially close to zero, as they vanish for Hermite polynomials. 
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Appendix C. Deformation equations for orthogonal polynomials 

In this Appendix, we prove the deformation equations (1.36) and (2.2). These equations are 

in fact quite general, and we present the proof for a general class of orthogonal polynomials. Let 
{Pk{^)}kLo be a system of monic polynomials satisfying the orthogonality condition 

/ pk{x)pj{x)e'''''^ djiix) = hkSjk, (C.l) 

JR 

where d(i{x) is any measure on M such that the system of orthogonal polynomials exists. We 
consider deformations of this system with respect to the parameter a. Let us write the three 
term recurrence equation, explicitly noting the dependence of each recurrence coefficient on the 
parameter a: 

xpk{x) = Pk+i{x) + Ak{a)pk{x) + Bk{a)pk-i{x) , Bk{a) = . (C.2) 

Notice that, since the polynomials pk are monic, -^pk{x) is a polynomial of degree strictly less than 

k, and thus its integral against pk{x)e~"-^ diJ,{x) is zero. Thus if we differentiate (C.l) with respect 
to a in the case j = k, apply the three term recurrence twice and integrate, we obtain 

h'k{a) = -hk{Bk+i + Ak + Bk), (C.3) 

or equivalently 

— log hk = -Ak , (C.4) 

da hk hk-i 

where we have suppressed the notation which explicitly indicates dependence on a. 

Let us use Ckj to denote the coefficient of the x^ term in the polynomial Pk{x), so that 

Pk{x) =x^ + Cfc,fc_ix"-i + Ck,k-2x''-^ + ■■■ . (C.5) 

These coefficients depend on the parameter a, and by matching the coefficients of the x^ term in 
(C.2), we see that 

Ak{a) = Ck,k-i - Cjfc+i,fe . (C.6) 
To arrive at a deformation equation for Ak consider (C.l) with j = k — 1. Differentiating with 
respect to a and disregarding the term for which the integral vanishes gives 

d 

-g^{Pk{x))pk-i{x) - x'^pk-i{x)pk{x) 
Applying the three term recursion twice and integrating, we obtain 

(^-Q^Ck,k-i^ hk-1 = Akhk + Ak-iBkhk-i ■ (C.8) 

Combining (C.6) with (C.8) both as it is written and with k t-^ k + 1, we find 

A'kia) = Bk{Ak + Ak-i) -Bk+i{Ak+i+Ak). (C.9) 
We now use (C.3) and (C.9) to integrate (C.3) once more, obtaining 

■^loghk = h+i- Ik, (CIO) 



e-"^'d/x(x) = 0. (C.7) 



where 



Ik = Bk(^Bk+i + Bk-i + {Ak-i + Ak)^^ . (C.ll) 
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It follows that the sum 

fc=0 

telescopes and its value is /„ — Iq. But /q = 0, and thus the sum (C.12) is simply /„. After a 
change of variable, this proves (1.36). 

We now prove (2.2). In the case that the measure of orthogonality is even, the recurrence 
coefficients vanish, and we have 

Ik = BkBk+i + BhBk^i (C.13) 

and 

log /l2fc = E B2k+lB2k+2 - B2kB2k-l ; E log /l2fc+l = ^ -B2A;+2-B2fc+3 " B2kB2k+l 

k=0 k=0 k=0 k=0 

(C.14) 

which are again telescoping sums, and we obtain (2.2) after a change of variables. 
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